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Abstract 

In unparticle physics, operators of the conformal sector have self-interactions, and these are 
unsuppressed for strong coupling. The 3-point interactions are completely determined by conformal 
symmetry, up to a constant. We do not know of any theoretical upper bounds on this constant. 
Imposing current experimental constraints, we find that these interactions mediate spectacular 
collider signals, such as pp — > U — > UU — ► 7777, 77ZZ, ZZZZ, ^yyl + U, ZZl + U, and 41, with 
cross sections of picobarns or larger at the Large Hadron Collider. Self-interactions may therefore 
provide the leading discovery prospects for unparticle physics. 

PACS numbers: 12.60.-i, 11.25.Hf, 14.80.-j, 13.85.-t 
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I. INTRODUCTION 



In unparticle physics, the standard model is extended by couplings to a conformal sector 
through interactions of the form 00<^m, where O is an operator of the conformal sector, and 
Osm is a standard model operator The conformal sector may be weakly-coupled [2J or 
strongly-coupled [3], but in all cases, the conformality of the new sector leads to effects that 
cannot be explained in terms of standard particle states. 

To date, all unparticle studies are based on two key elements: the unparticle phase space, 
and the unparticle propagator. Conformal invariance fixes unparticle phase space [l[ , which 
enters processes with unparticles in the final state, such as ff — > fU. Similarly, conformal 
invariance dictates the form of the unparticle propagator jj, [^], which determines virtual 
unparticle contributions to processes such as // — > U — > //. The forms of the unparticle 
phase space and propagator imply that unparticles do not behave as standard particles, but 
are more aptly interpreted as fractional numbers of massless particles [l[ or as collections of 
particles with a particular distribution of masses 0] . These results are valid if the couplings 
of unparticles to the standard model are all non-renormalizable. If unparticle couplings 
include the super-renormalizable operator Oh 2 , where h is the standard model Higgs bo- 
son, electroweak symmetry breaking breaks conformal invariance 0, [1] ■ This modifies the 
unparticle propagator and implies that unparticle physics may be probed only in a narrow 
conformal window, typically at energies between 10 GeV and 1 TeV 0, @|. Such energies 
are best probed at high energy colliders, and many studies have investigated the collider 
implications of unparticles p], 13, 0, @, [Hi Ell, E3] ■ Models with conformal breaking may also 



share characteristics with hidden valley models [13j, [14] . 

Here we study a qualitatively new effect: unparticle self-interactions. Such 3- and higher- 
point interactions are always present in conformal theories, and mediate processes such as 
gg — > U — > U ■ ■ -U, with two or more unparticles in the final state. In the most interesting 
cases with strongly-coupled conformal sectors, the creation of additional high unparticles 
in the final state does not suppress the rate. Multi-W production therefore differs from all 
known examples, such as gg — > g — > g ■ ■ ■ g and gg — > 7 • • • 7, where the rate is reduced with 
the addition of every high particle. 

In this paper, we focus on 3-point self-interactions. These are the natural starting point 
for several reasons. First, unlike 4- and higher-point self-interactions, 3-point interactions 
are completely constrained by conformal invariance, up to a constant. Second, 3-point 
interactions are the leading order at which induced signals may be nearly background free. 
For example, although the IA propagator induces signals like pp — > hi — > 77, 3-point U self- 
interactions mediate pp —>■ U —>■ UU — > 7777, 77ZZ, ZZZZ, / -f-fl + l~, ZZl + l~ ,41, and many 
other spectacular signals through subprocesses such as the one shown in Fig. [TJ For this 
reason, the most promising signals for unparticle discovery at colliders may in fact be those 
that are induced by 3-point unparticle interactions. (Note that, although the production of 
3 or more unparticles may also be unsuppressed, the requirement that they convert back to 
standard model particles to be visible does, in fact, imply that these are sub- dominant.) 

In the following sections, we start from the 3-point correlation function in position space 
and convert it to a form useful for Feynman diagram calculations. We then calculate collider 
rates, taking the 47 signal as our example. Our final result is determined up to the constant 
entering the 3-point correlation function. As far as we know, there is no theoretical upper 
bound on this constant — it is bounded only by current experimental constraints. Applying 
these constraints, we find that the prospects for signals at the Large Hadron Collider (LHC) 



2 



Pi 




FIG. 1: The process gg — > 7777 mediated by unparticle self-interactions. 

are truly spectacular - - the four-body final states given above may have cross sections 
of picobarns or larger, leading to obvious signals visible in the first year of the LHC We 
also determine the predicted kinematic distributions from unparticles, which are insensitive 
to overall rates. These distributions provide testable predictions that could be used to 
distinguish multi-unparticle production from other possibilities for new physics. 



II. 3-POINT CORRELATION FUNCTION 



We begin by assuming that the standard model is coupled to the unparticle sector through 
a scalar operator O. For a scalar, unitarity requires dimension d >__1, but there is no upper 



bound [15|, [16|. Motivated by prominent supersymmetric examples [3j, we consider the range 
1 < d < 2. Modifications are required for vector and tensor operators 15|, [l6|, [17|, but we 
expect that our primary conclusions apply to these cases as well. 

Conformal invariance constrains the 3-point O correlation function to be 
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where C' d is a constant determined by the unparticle self-coupling strength. In momentum 
space, it is then 
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C d is 



where the last equality makes use of the unparticle 2-point correlation function 
determined in terms of C' d ; we choose to express our results in terms of C d - 

The 3-point correlation function is the product of 3 propagators from to x, x to y, and 
y to 0. It is therefore not surprising that, in momentum space, it takes the form of a loop 
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integral. We may therefore make use of the standard techniques available for simplifying 
loop calculations. Using Feynman parameters for Eq. ([2]), the 3-point correlation function 
becomes 



(0\O( Pl )O(p 2 )O\ Pl +p 2 )\0) = -*(-!)» C„ (iy 2 FJd t d\ , ( 3 ) 



where 
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n = 6 — y and A' = xia^Pi/s + ^la^P^/ 5 + ^2^3, with s = (pi + V2) 2 ■ Fy(pl/ s iP 2 / s ) is 
plotted in Fig. H 
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III. BOUNDS ON UNPARTICLE INTERACTIONS 



To investigate the phenomenological implications of the 3-point correlation function dis- 
cussed in Sec. [TTJ we must introduce unparticle couplings to standard model particles. Here 
we present here the relevant couplings, and determine the existing constraints on them and 
the prospects for probing them at the LHC. The results here are independent of unparticle 
self-interactions, but are of interest in their own right. In addition, in Sec. II VI we will present 
results for a reference value of A4, which we choose based on the results we derive here. 

A scalar unparticle may couple to two gauge bosons or two fermions through the La- 
grangian interaction terms 



(5) 



leading to the Feynman rules 
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where A 4 is some high scale characterizing these non-renormalizable interactions, v is the 
Higgs vacuum expectation value, e is the proton charge, and c 7 , c g , and c{ are constants. 

Following convention, we choose an O normalization by specifying the unparticle propa- 
gator 0, Q 



scalar unparticle propagator : iBdO(q°) 0(q 2 ) (q 2 — yU 2v '' 1 



(7) 



where 
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Here the modified propagator suggested in Ref. [?J is used to take into account the break- 
ing of conformal invariance at a scale fi by unparticle couplings to the Higgs boson. 
A more detailed analysis was performed in Ref. [8] by considering a deconstructed ver- 
sion of the unparticle-Higgs coupling. That approach led to a propagator of the form 

(q 2 ) 2 ~ d — v 2 ((J>u) 2 ~ d / (q 2 — m 2 ) , where v and m h are the Higgs vacuum expectation value 
and mass, respectively, and \iu is a scale related to the Higgs-unparticle coupling. 

We will be most interested in unparticle couplings to two photons. This coupling, in 
conjunction with unparticle couplings to gluons and quarks, mediates processes gg, qq — > 
U — > 77, which have been studied previously in Ref. With our coupling conventions, 
these processes have differential cross sections 
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FIG. 3: Bounds on the unparticle scale A4 from pp — > IA — > 77 using existing data from the Teva- 
tron (solid), and the projected 5a discovery reach at the LHC for integrated luminosity 100 fb _1 
(dotted). For comparison, existing bounds from e + e~ — > U — > fj> + fi~ @] (long dashed) and the 
unitarity bound on WW scattering lCj] (short dashed) are also shown. See text for details. 



where we have averaged and summed over initial and final state colors and polarizations, 
but have not yet included the factor 1/2 to account for the two identical photons in the final 
state. 

These processes are bounded by existing data, most stringently by the diphoton invariant 
mass distribution from the CDF Collaboration at the Tevatron [l8( , based on an integrated 
luminosity of 1.2 fb _1 at v/i = 1.96 TeV. The events in this distribution have two central 

12 m 

photons with \rj\ < 1.04 and transverse momenta pj, > 15 GeV. To obtain a lower bound 
on A4, we simulate the signal by adopting the following procedure here and in all analyses 
described below. We set c g = c 7 = 1 and e 2 (c{) 2 = 2ir, following the convention justified 
in Ref. use the CTEQ5L parton distribution functions with a factorization scale \l$ 



\ph 19], and evaluate the cross section with the Monte Carlo program VEGAS 20] . We 
then impose the identical rj and pt cuts given above, and further require m 77 > M t h, where 
m 77 is the diphoton invariant mass, and M t h is a threshold mass chosen to maximize the 
sensitivity to the signal. At the Tevatron, unparticle production through qq dominates over 
gg, and the most stringent bounds are typically achieved for M t h ~ 350 GeV. The 95% CL 
lower bound on A 4 is then derived following Ref. [2l|; these results are given in Fig. For 
comparison, in Fig. [3] we also present bounds on A 4 given previously in the literature from 
e + e~ — » IA — > from LEP/SLC [9] and the unitarity bound on WW scattering Tol . 

We have also determined the discovery reach for pp — > U —>■ 77 at the LHC, assuming 



100 fb 1 of data. We evaluate the prompt diphoton background using PYTHIA 6.4 22 



impose kinematic cuts \rj\ < 2.5 and p^ > 25 GeV, and again optimize by varying the 
threshold energy M t h- We estimate a 5a discovery reach by requiring S/ \[B > 5 and S > 5, 
where S and B are the number of signal and background events. This discovery reach is 
also presented in Fig. [31 

These results have omitted important effects. Notably, in determining the LHC reach, 
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we have omitted the background from misidentification of jets and electrons. Nevertheless, 
we may conclude that present limits from the Tevatron and LEP/SLC are comparable, and 
imply lower bounds on A4 of ~ 2 TeV for d = 1.1, dropping to below a TeV for d = 1.9. 
We will adopt the reference value A 4 = 1 TeV below. At the LHC, the sensitivity is much 
higher, with the 77 signal probing A 4 in the multi-TeV range. Note that here, and in the 
rest of this work, we have taken the conformal breaking scale to be /i = in deriving our 
results. However, we have found that our results are insensitive to this choice, provided 
\i < M t h. In general, the modifications have the property that they are small for large 
q 2 . For example, for q 2 3> /x 2 , the propagator we used (Eq. (j7])) is of the form (q 2 ) d ~ 2 with 
corrections of order fi 2 /q 2 . In this limit the propagator obtained from Ref. [8j] has corrections 
of order (/i 2 / q 2 ) 2 ~ d v 2 / (q 2 — ml). We will assume henceforth that the breaking of conformal 
symmetry is small in the kinematic region we work in so that the correction terms can be 
ignored. 



IV. FOUR PHOTON EVENTS 

We now turn to new processes mediated by unparticle self-interactions. As noted in Sec. [Ij 
the 3-point function is particularly interesting, and, depending on the unparticle interactions 
with the standard model, may mediate a variety of processes leading to spectacular 4-body 
final states. As an example, in this section we consider the four photon signal shown in 
Fig. EJ 

The cross sections for gg,qq — > U — > UU — > 7777 are completely specified, given the 
Feynman rules of Eqs. ([3]) and ([6]). With momenta as labeled in Fig. [TJ we find 
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where we have averaged and summed over initial and final state colors and polarizations, but 
have not yet included the factor 1/4! to account for the four identical photons in the final 
state. To evaluate the parton- level cross sections, we use the results of Sec. [Til for the 3-point 
correlation function. Setting c g = c 7 = 1 and e 2 (c{) 2 = 2n as above, and integrating over 4- 
body phase space with the Monte Carlo program VEGAS [20(] , we find that the parton- level 
cross sections may be written as 
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where the dimensionless proportionality factors /| and /J are given in Table [TJ 
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3.6 x 10~ 7 
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LHC maximum cross section [fb] 


4300 


9600 


1.8 x 10 5 


8.4 x 10 6 


ILC (# e __^ [fb] 


1.0 x 10~ 6 


4.7 x 10~ 7 


6.0 x 10~ 8 


8.0 x 10" 9 


ILC maximum cross section [fb] 


160 


250 


810 


1900 



TABLE I: Dimensionless parton-level proportionality factors f^' q , and unparticle 47 reference cross 
sections at the Tevatron (y/s = 1.96 TeV), LHC (y/s = 14 TeV) and ILC (y/s = 1 TeV). These 
reference cross sections assume A 4 = 1 TeV and Cd = 1; the actual cross sections scale with 
C^/A d d . Upper bounds on this combination of parameters from existing data at the Tevatron are 
also listed, as are the resulting maximum cross sections possible at the LHC and ILC. 

A. Bounds from Tevatron 

The 47 unparticle signal is bounded by searches at the Tevatron. The DO collaboration 
has searched for the inclusive production of multi-photon final states j23[. Events are selected 
with three or more photons in the central calorimeter (|^| < 1.1) and ^-ordered cuts on 
their transverse energies: E^ 2 ' 3 ^ > 30, 20, 15 (, 15) GeV. The dominant backgrounds are 
diphoton production with additional initial state radiation (ISR) photons and events in 
which jets or electrons are misidentified as photons. No excess of events above the standard 
model prediction was found in the Tevatron data with integrated luminosity 0.83 ±0.05 fb -1 
collected during 2002-2005. 

This search result thus sets upper bounds on combinations of the unparticle self- 
interaction strength Cd and the energy scale A 4 . To derive these bounds, we again set 
c g = c 7 = 1 and e 2 (c{) 2 = 2n. As with the 77 signal discussed in Sec. IIII[ at the Tevatron, 
the qq contribution dominates over that from gg, since pp collisions provide a large density 
of anti-quarks, and at y/s = 1.96 TeV, large Bjorken x is required. 

We impose the kinematic cuts on rapidities and transverse energy as given above. Esti- 
mating the background event rate to be C cr^ 4 m C (J2 q Q 2 ) («/tt) 2 a qqL2-y ~ 0-83 fb -1 x 
(0.002) 2 x 0.66 pb ~ C(10- 3 ), we find that it is negligible. The 95% CL upper limit of 3.04 
events for zero background and zero events observed therefore becomes the 95% CL bound 

r ,( 1 V% 3.04 
d yA 4 [TeV] J ~ 0.83 ftr 1 ' K } 

where = cXg| f _^ 47 + cr^ f _^ 47 is the total reference cross section, determined by setting 
A 4 = 1 TeV and Cd = 1. The resulting upper bounds on C^/A 3 ^ are given in Table [fl 
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FIG. 4: Maximum possible cross sections for 4j production through unparticle self-interactions, 
given existing bounds on 47 events from the DO Collaboration at the Tevatron 23]. 



B. Prospects for LHC and ILC 

We now determine the LHC cross sections for 47 production mediated by unparticle self- 
interactions. Using VEGAS j20|, we Monte Carlo simulate 4 to 10 million events for several 
values of d. We require the photons to have rapidity \q\ < 2.5 and transverse energies 
£,1,2,3,4 ^ ^ Q e v. The resulting reference cross sections for A4 = 1 TeV and 

Cd = 1 from gg and qq initial states are given in Table [B 

Of course, the actual cross sections scale with C\. We do not know of any way to bound 
Cd theoretically, or, indeed, to specify a "typical" value for Cd- Lacking theoretical guidance, 
we simply impose existing experimental bounds. Although bounds on unparticles have been 
discussed at length in the literature, these constraints, including those derived and presented 
in Sec. 1111} do not bound Cd- For a direct bound, we must therefore turn to the Tevatron 
results derived above, which, of course, constrain Cd and A 4 in the combination identical to 
that which enters the LHC cross section. Assuming these parameters saturate the bound of 
Eq. f|T5|) . we find the maximal cross sections presented in Table |T]and plotted in Fig.Hl These 
cross sections are extraordinarily large, ranging from picobarns for d — 1.1 to nanobarns 
for d = 1.9. If the unparticle self-interaction is anywhere near the largest values allowed by 
current experimental constraints, this spectacular signal will be discovered very early at the 
LHC. 

The px distributions of the leading, next-to-leading, next-next-to-leading, and fourth 
photons in unparticle 47 events are given in Fig. We find that the px spectra, even for 
the 3rd and 4th photon, are remarkably hard. Given this, the standard model background 
of diphoton production with two ISR photons can be eliminated with pt cuts with little 
effect on the signal. The dominant backgrounds will be from events with misidentified jets 
or electrons, but even these may be reduced significantly with hard pt cuts without large 
reduction in the signal. It would be very interesting to determine the extent to which this is 
validated by a realistic detector study. Last, we note that the shapes of the px distributions, 
along with other kinematic information, are, of course, independent of Cd and A 4 . They may 
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FIG. 5: The pt distributions for p^-ordered photons in unparticle 47 events (from the dominant 
gg initial state only) at the LHC for d = 1.1 (top) and 1.9 (bottom). 

therefore be used to identify unparticles as the source of these signals, and to distinguish 
them from anomalies predicted by other frameworks for new physics. 

If the coupling of scalar unparticles to electrons is of the same size as to quarks, 47 event 
rates at the International Linear Collider (ILC) can be calculated with Eq. ffl~4j) multiplied 
by 3. ILC reference cross sections and maximal possible cross sections are presented in 
Table fl] and Fig. |U Interestingly, the 47 cross section at a 1 TeV ILC is about two orders 
of magnitude larger than at the Tevatron. This is because the parton densities become less 
than one for x > 0.3, so the Tevatron pp — > 47 cross section is dominated by the contribution 
from VI ~ 300 GeV. 



V. DISCUSSION 

In this work, we have investigated a new feature of unparticle physics, namely, the self- 
interactions of unparticles. Such interactions are necessarily present in conformal theories, 
are unsuppressed for strongly coupled conformal sectors, and introduce a large range of new 
phenomena not studied previously. 

As an example, we have investigated the 3-point correlation function in detail. This is 
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completely specified by conformal invariance, 1 up to a constant Cd, and mediates a host 
of processes, such as pp — »■ U — > UU — > 7777, 77ZZ, ZZZZ, , yyl + l~ , ZZl + l~ ,41 at the 
LHC. We know of no upper bound on Cd- We have therefore imposed only the existing 
Tevatron bounds, and have found that, for the example 47 signal discussed here, the allowed 
cross sections at the LHC are enormous, ranging from picobarns for unparticle dimension 
d = 1.1 to nanobarns for d — 1.9. Such signals could emerge very early after the LHC turns 
on, and imply that the 4-body final states could be by far the most promising modes for 
unparticle discovery at colliders. Given these observations, it would be very interesting to 
perform similar analyses for the other final states, as well as to carry out realistic detector 
simulations to determine the extent to which jet and lepton misidentification masks these a 
priori spectacular signals. 

On the theoretical side, the most pressing issue is to determine what values of Cd are 
possible or natural. Conformal invariance by itself does not determine the constant, and 
we are unaware of any other consistency conditions that could bound Cd- It would be 
very interesting to see if the AdS/CFT correspondence can shed light on this issue. In the 
absence of theoretical bounds, Cd can be large, and pp — ► UU — ► 47 and related 4-body final 
states may have cross sections larger than pp — > U -^77 and related 2-body states. Such 
a conclusion is highly counter-intuitive, but perhaps possible, given the other surprising 
properties of unparticles discovered so far. 
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